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1. Introduction 


The problem of nonlinear Helmholtz resonators has been addressed by 
many previous investigators. The character of the nonlinearity and various 
additional aspects of resonant and nearly resonant oscillations have been 
discussed by Ingard [1], Bies and Wilson [2], Gill and Bhaduri [3], Barthel 
[4], Borisova et al. [5] and Ochmann [6]. However, there still doesn't seem 
to exist a systematic investigation of the relative importance of a variety of 
effects, including those of thermoacoustic Stokes boundary layers at the 
walls of the throats, a non-zero mean flow through the resonator, the 
conductance of the ends of the throats and higher-order frequency correc- 
tions to the *Helmholtz approximation". 

In the limit of the Helmholtz approximation spatial phase variations 
within the individual elements (throats.and volumes) of a resonator become 
negligibly small. In this limit a resonant oscillation may be regarded as a 
periodic exchange process that involves conversion of potential energy from 
pressurized volume elements into kinetic energy of moving gas columns in 
throat elements, and vice versa. Due to radiation of sound, boundary layer 
dissipation and internal jet dissipation a (small) fraction of the acoustic 
energy, stored in a resonator, gets lost per cycle and must be compensated 
by an excitation mechanism to sustain an equilibrium oscillation. Typically 
such an excitation is caused by a pressure oscillation at an outer end of a 
throat element. To obtain an accurate prediction of oscillation amplitudes 
all the aforementioned mechanisms of excitation and attenuation must 
usually be taken into account. Furthermore, lowest-order theories of nonlin- 
ear acoustic oscillations typically predict the resonance frequencies to 
first-order accuracy only. Higher-order frequency corrections are usually 
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Figure 1.1 
Arrangement of a Helmholtz resonator with two throats and flushing. 


taken into account, because most technical applications require knowledge 
of resonance frequencies that is of better accuracy than first order. 

For some technical applications a principal additional difficulty arises 
due to the necessity of permanently flushing the resonator, in order to avoid 
an overheating produced by the radiation of hot gas. The idea is to develop 
a design method which leads to reliable predictions of the resonance 
frequency, the width of the resonance band and the attenuation power of a 
resonator element. Although, the width of the resonance band can be 
changed by changing the geometry of the Helmholtz resonator, it is clear 
that the relative width of the resonance band is always fairly small. For this 
reason Helmholtz resonators are very effective attenuators for narrow 
frequency intervals. However, they are not at all suitable for broad band 
sound attenuation. For a variety of configurations and conditions the 
properties of Helmholtz resonators have been investigated with the help of 
the experimental results reported in this paper. It will be seen that the effects 
of throughflow, the geometry of the ends of throats, higher order correc- 
tions of the Helmholtz approximation and Stokes boundary layers in the 
throats are usually of fundamental importance for an accurate theoretical 
prediction of the properties of a Helmholtz resonator. 

The situation presently considered is depicted in Fig. 1.1. 


2. Resonance frequency in the limit of negligibly small boundary layer friction 


The Helmholtz resonator shown in Fig. 1.1 consists of the three tube 
elements designated by 7, II and MI. Assuming that the diameters and 
lengths of the tube elements are very small compared to the relevant 
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acoustic wave lengths and that the attenuation of acoustic waves within the 
tube elements due to boundary layer friction is negligibly small, the down- 
stream and upstream running acoustic waves (i.e. Riemann invariants) are 
defined by the expresssions 

pe Ea plc (2.1) 

Qc Qc 

where Au and Ap denote velocity and pressure disturbance, o refers to the 
density and c to the sound speed. All quantities which are constant within 
the tube elements of the Helmholtz resonator are marked with suffices 7", 
"IT", " IIT" respectively. Hence, the mean flow Mach number is denoted by 
M,, for example. The discussion in this section is restricted to the small 
amplitude limit of acoustic oscillations. Hence, linear acoustic theory is the 
adequate tool for the subsequent consideration. The reflection condition at 
the inlet of the first tube of the resonator can be defined by 


_1=GM, , 
1-0M;TU 
where “1” refers to the location in the resonator (see Fig. 1.1) and £ is the 
loss coefficient of the tube inlet; £, = | would correspond to a loss free inlet 
with a static pressure drop of one dynamic head across the tube inlet. 


As only the relative amplitudes of the acoustic waves are relevant for the 
present linearized analysis, we introduce the normalization 


fh 21-6M; and g)=1+¢,M;. (2.3) 


Harmonic Riemann invariants f and g at the locations 1 to 6 are related by 
the expressions 


fi (2.2) 
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where A, L and o denote the cross-sectional area of a tube element, its 
length and the angular frequency of the harmonic waves respectively. 

While the phase shifts defined by the relations (2.4), (2.6) and (2.8) are 
obviously related to the lengths of the tube elements and the mean flow 
speed M : c, the relation (2.5) is obtained on the basis of the assumptions 
that the static pressures at the locations 2 and 3 are equal (fully separated 
free jet) and that the phase shift between these two locations is negligibly 
small. The latter assumption also holds for the relation (2.7). However, the 
static pressure drop between the locations 4 and 5 is taken into account. For 
a loss-free flow between these two locations the coefficient (jjj would assume 
the value 


An V 
m=1—-|—~—}. 2.9 
( a) (23) 


Lec 0 
| | M (1+ Mg): c | j| | (2.8) 


L 


To obtain the pressure oscillation at location 6 we evaluate the term 
$( fe — gg). Within the “Helmholtz approximation", which is discussed sub- 
sequently, the Mach number corrections contained in the phase exponents 
of (2.4), (2.6) and (2.8) can be neglected. Combining the relations (2.3) to 
(2.8) and neglecting the Mach number corrections in (2.4), (2.6) and (2.8) 
leads to the expression (A1) given in the Appendix. 

Now we restrict our consideration to the lowest resonant frequency and 
make use of the “Helmholtz approximation". In this case the expression 
(Al) can be simplified by introducing the following hierarchy of orders: 


e) =o) tne 
C c C 


(2.10) 


and 


O(M7) = O(M m) S 8. (2.11) 
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For a given amplitude of the pressure oscillation at location 6 a necessary 
condition for resonance is that the amplitude of the pressure oscillation in 
the vessel /7 reaches a local maximum. However, because we have already 
fixed the amplitudes of the waves by the normalization (2.3) we have to 
keep the amplitude of the pressure oscillation at location 6 variable. For this 
reason the corresponding necessary condition for resonance is that the ampli- 
tude of the pressure oscillation at location 6 reaches a local minimum. 
Making use of the substitution 


1 
5 (Je — 8g.) =E+M, where € and 7 are real, 


this condition can be written in the form 
0 2 2 
—(E’+ =(. 2x 
5 (C q^) (2.12) 


Inserting the expression (Al) into the condition (2.12) and making use of 
the hierarchy of orders as defined by (2.10) and (2.11) leads to lowest order 
to (A2) as given in the Appendix. 

In the limit M;, Mjj;-—0, the resonance condition reduces to the 
condition that the underlined term within the first square bracket of (A2) 
vanishes. Thus we have 

wo? 1 [Ar Am 

eo 4 [2 T Lin 
which represents the well-known Helmholtz formula for a resonator with 
two throats. 

From a practical point of view the following special case is particularly 


important 


Ar | 
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It is quite remarkable that in this case the general condition (A2) reduces to 


i where V = Lg 3 An, (2.13) 


2 
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eU FEES! 
which is the simple Helmholtz formula for a resonator with one throat only. 

In the general case where no further simplifications are justified the 
resonance frequency is obtained by solving equation (A2) for (c/c). 


(2.15) 


3. The attenuation of waves due to the thermoacoustic boundary layer 


Assuming that the mean velocity in the tube elements (i.e. the throats 7 
and JII) is very small and that the thicknesses of the boundary layers are 
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small compared to the tube radii, leads to an exponential attenuation of 
harmonic waves: 


f=a- exp i — — — 3 
and (3.1) 
geb ex (io T = + ax), 
c 


where |a| « œ/c. 
Following Rott and Hartunian [7] (see also Chester [8] and Keller [9]) 
we find 


| (3.2) 


where v denotes the kinematic viscosity of the air, D the diameter of the 
tube, y the ratio of specific heats and Pr the Prandtl number. 

Hence, the attenuation factor of a harmonic wave running from one end 
of a throat of length L and cross-sectional area A to the other is 


exp[ — Re(x) : L] = exp} - E PE = [ zl (3.3) 


where Re(«) refers to the real part of a. 

As long as the mean flow speed in a tube does not exceed the order of 
magnitude of the amplitude of the velocity oscillation, the viscous effects 
due to the oscillatory part of the velocity are not affected to first approxima- 
tion by the mean flow speed. 

Hence, we can restrict our discussion to the case of zero mean flow 
speed. To fix ideas we consider a pair of upstream and downstream running 
waves of the form (3.1) in a tube of length L and diameter D. Furthermore, 
we set the origin of the coordinate system such that x = —iL and x = +4L 
correspond to the two ends of the tube (see Fig. 3.1). 

For reasons of symmetry we can now set b — a. 

Combining (2.1) and (3.1) yields 


Au = ; (f+g)=a-e™. cosh( 12 + ax) (3.4) 


— => (f8) = -a e. sion" + 3! (3.5) 


Vol. 46, 1995 The use of Helmholtz resonators as sound attenuators 303 


© 


Figure 3.1 
Coordinate system in a tube element. 


Making use of the Helmholtz approximation as defined by (2.10), the 
expressions (3.4) and (3.5) can be simplified: 


A ji 
Au =a: e, La —aq-:e^t. E + ax | (3.6) 
Qc c 
Hence, we obtain (see locations marked in Fig. 3.1) 
; iL 
Ap, — Ap —a : e! - [SE esL | oc (3.7) 
Introducing now the boundary layer parameter 
y—1 
DA 
Pr!? 2v 
- 3.8 
s D 7 (3.8) 
the equations (3.2), (3.6) and (3.7) yield 
ö 
Ap, —Ap,—-[l-cs]:o:- L- z ^u +s:0: L-: c Au. (3.9) 


Ot E 


The first term on the right hand side of equation (3.9) obviously corre- 
sponds to the inertia of the gas column in the tube, including a correction 
due to the boundary layer, and it is non-dissipative, whereas the second 
term accounts for the viscous effect on the pressure difference along the 
tube. 


4. The Helmholtz approximation of the equations of motion 


We return now to consider the complete resonator as shown in Fig. 1.1. 
Within the Helmholtz approximation, which was already discussed at the 
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end of Section 2, the condition for mass conservation can be expressed as 
follows: 


QV = (U; + ur) -Q ` Aj (Ur + ug) ` 0 ^ Agr, (4.1) 


where U; and Uy denote the mean values of the velocities in the tube 
elements 7 and 7/1, respectively, and u, and uyy refer to the departures of the 
velocities from their mean values. Within the same approximation (4.1) 
implies that 


U;A, ra OA =0 (4.2) 
and 
QV —uj: g^ Aj —umr ` Any. (4.3) 


Denoting the pressure in the pulsating system near the end of the tube 
element JI] by p, and the pressure in the resonator vessel // by p, and 
ignoring in a first step the boundary layer effects, the equation of motion of 
the gas column in the tube element MI can be written in the form 


Ai Ly Qin = Amp, — Am ` Snr’ 2: (Ur + trm) ` JU sr F un) Cue Aug Pes 
(4.4) 


where a dot above a variable refers to a time derivative. 

The term on the left-hand side of equation (4.4) represents the force due 
to the acceleration of the gas column, the first and third terms on the 
right-hand side account for the force due to the different pressures in the 
resonator vessel and the pulsating system, and the second term on the 
right-hand side of equation (4.4) accounts for the force due to the accelera- 
tion of fluid from the resonator vessel into the tube element J or from the 
pulsating system into the tube element JI, depending on the sign of the 
velocity (U; + Up). In the case of loss-free acceleration ¢,,, is defined by 
(2.9) for positive velocities and £;; « 1 for negative velocities. For simplicity 
we assume that even for negative velocities hot fluid does not penetrate from 
the pulsating system into the tube element ZI. Hence, the density o does not 
change within the approximation considered. Here we assume that Cry 
assumes the same value for positive and negative velocities. However, a 
more detailed investigation shows that the general form of the second term 
on the right-hand side of equation (4.4) remains unchanged even for a 
strong dependence of (,, on the sign of the velocity, provided a loss 
coefficient is defined which corresponds to an average over an oscillation 
period. It should be pointed out that such an average strongly depends on 
the relative magnitudes of U, and uy. In situations of practical interest, in 
which the Helmholtz resonator is used as an attenuator, the aim is to keep 
(jy as small as possible. In this case “y; would usually not depend on the 
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sign of the velocity (U; + uz), because the two ends of the tube element JI 
would typically exhibit the same geometry. 

On the basis of the Helmholtz approximation (see (2.10) and (2.11)) it 
is easy to show that 


O(Quj) < (oüz) (4.5) 


and 


O(Q [Un + uml) < Oleum). (4.6) 


Differentiating (4.4) with respect to time and making use of (4.5) and (4.6) 
yields 


Lin’ ` ling = Pe — Gi Um |Unr + un — pe. (4.7) 


It should be pointed out that a mean-flow pressure drop in the tube element, 
which is produced by boundary layer friction, can be included by means of 
a suitable correction of “r. However, the effects produced by a Stokes 
boundary layer are different. A more general version of the equation of 
motion, which does account for the Stokes boundary layer, is obtained by 
including the corresponding terms according to (3.9) in equation (4.7): 


[1 + sud]: Lug Qr Up; = Po — Sar Q © Umr’ Um + ug 
— Sm ` Lig Q7 OO ^ gy — De (4.8) 


where 


y—1 ado TY 
su =| init \ 2:0 Any (4.9) 


For tube elements with non-circular cross-sections the diameter ‚/4A /r in 
(4.9) should be replaced by the hydraulic diameter. 

For the subsequent discussion we assume that the time-dependent part 
of the pressure in the pulsating system (near the end of the tube element 777) 
is a harmonic oscillation, i.e. 


Pe = De? <p i sin(œ ; [t p tol), (4.10) 


where “< »" refers to a time average, f is an oscillation amplitude and ¢, 
accounts for a shift of the origin of the time coordinate which will be chosen 
suitably below. Furthermore, we assume that the pressure fluctuation in the 
chamber upstream of the resonator is negligibly small. Inserting (4.10) in 
(4.8) yields 


[1 + Sy) © Lg: Qo Üm = Py — norme lU; T ur — Sm Lin 


Qr Og +@-p-cos(@ - [t — tj). (4.11) 
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A corresponding consideration leads to the equation of motion for the tube 
element T: 


[+s]: Leo ü= — tr eo úr |Ur+u|—-s7 Leo ùr p. 
(4.12) 


Within the present approximation it is justified to use the following relation 
(which is strictly valid for isentropic flow): 


p. = 076. (4.13) 
Combining (4.3) and (4.13) leads to 


p= Aruroc? = Animo’ 
: V V ` 


(4.14) 


As will be shown subsequently, a harmonic pressure oscillation of the 
pulsating system also leads to harmonic velocity oscillations in the tube 
elements and a harmonic pressure oscillation in the resonator vessel within 
the approximation considered. 

The time shift tọ appearing in the expression (4.10) is now implicitly 
determined by choosing the following expression for the velocity oscillation 
in the tube element ZT: 


Urg = Ür t sin(@ ` D). (4.15) 
For the velocity oscillation in the tube element J we write 
ur = —üj:sin(o [t — x]). (4.16) 


The symbol “^” again refers to amplitudes, and the phase shift associated 
with the constant t will have to be determined subsequently. 

Inserting (3.14) in (4.11) and (4.12), in order to eliminate p,, leads to 
Am’? Aroco, 


(i + Sr) ` Lg: Q` Üm + y "Un V Ur 


= —óm'Q'Up Un; + Un] — Sm ^ Lin Q^ O ^ üm 
+o : p-cos(@ : [t — tol) (4.17) 
and 


Arge _ Anr E 


(1 +5) Le ü+ V Ur y ` Ur 


= —¢r 0 ` úr: |Ur + ur — sr’ Lrg o ty. (4.18) 
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5. Properties of the oscillation equations at and near resonance 


With the equations (4.17) and (4.18) we have now obtained a complete 
system to describe small amplitude nonlinear oscillations of the resonator 
shown in Fig. 1.1 at and near resonance. 

If the right-hand sides of these equations were equal to zero, the system 
would become linear, and it would describe a neutrally stable coupled 
oscilation. The third term on the right-hand side of equation (4.17) 
accounts for the excitation, the first terms on the right-hand sides of the 
equations (4.17) and (4.18) describe the nonlinear attenuation due to free jet 
dissipation and the second terms account for the linear attenuation due to 
the Stokes boundary layer. 


Special case: the simple Helmholtz resonator 


In order to simplify the oscillation equations we first need to study 
certain basic properties of the equations. For this purpose we do not need 
to consider the relatively complicated general case. It is sufficient to investi- 
gate the envisaged properties in the case of the simple Helmholtz resonator 
as illustrated in Fig. 5.1. 

Setting A; — 0, to — 0, U; = Un — 0 and dropping the suffix “IIT”, the 
system ((4.17), (4.18)) is reduced to a single equation: 

cd ud 
V 
-—(:ocücl|u]—s: L: o: cou o p coslo : f). (5.1) 


This equation describes the oscillation of an ordinary Helmholtz resonator 
without throughflow. For a harmonic excitation as described by the third 
term on the right-hand side of (5.1), the solution u can be written as a 


Li Lu 


All Ap, =-Pp-sin (6 -t) 


Figure 5.1 
Helmholtz resonator with one throat only (closed resonator). 
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Fourier series: 
u= Y, [a,: cos(ot) + b, : sin(ot)], (5.2) 
n-l 
where a, and b, are constants. If the right-hand side of (5.1) were equal to 
zero the oscillation would again become neutrally stable. The corresponding 
resonant frequency is 


A 1/2 


Hence, at resonance (i.e. for c = ox) 


A 2 
(149) Lt Tu =0, (5.4) 
where 
u,-—d,:cos(o : t) +b, sin(o - t). (5.5) 


For this reason the third term on the right-hand side of equation (5.1) has 
to be compensated by the first and second terms on the right-hand side of 
(5.1) alone, if @ = œw. Thus, if the effects due to the Stokes boundary layer 
are very small, we have 


p 1/2 
max(O(a;), O(b;)) = (4) (5.6) 
On the other hand, if the attenuation due to the Stokes boundary layer is 
much stronger than the nonlinear attenuation, we obtain 


p 
max(O(a,), O(b,)) = (. TIE =} (5.7) 
In the latter case it is obvious that the higher harmonics in (5.2) are very 
weak compared to the fundamental harmonic. However, it can also be 
shown that the same is true for weak effects of the Stokes boundary layer: 
Inserting (5.2) in (5.1) and ignoring the effects due to the Stokes boundary 
layer (i.e. s x 0) we obtain 


s 

max(O(a,), Ob) = | p (5.8) 
Hence, 

max(O(a;), O(b;)) < max(O(a;), O(b,)), (5.9) 
if 

H x ; (5.10) 
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In a situation of practical interest the left-hand side of (5.10) would 
typically be two orders of magnitude smaller than the right-hand side. 
Similarly, it can been shown that the complete departure of u from its 
fundamental harmonic remains small if (5.10) applies, i.e. 


o( Y. [an coslo - t) + b, ` sin(a - at) < Olu). (5.11) 
n=2 

As the amplitude of u decreases rapidly as jw — wo| increases, and the first 
term on the right-hand side of (5.1) becomes very small compared to the 
left-hand side of (5.1), the relation (5.11) remains valid for œ # wo. 

To summarize the discussion we may conclude that the solution u of 
equation (5.1) is nearly sinusoidal at and near resonance, provided the 
amplitudes of the pressure oscillations remain small according to the 
criterion (5.10). 


General case 


For rather obvious reasons this general statement remains to be true if 
we extend the discussion again from the oscillation in an ordinary 
Helmholtz resonator as described by equation (5.1) to the more general case 
discussed in Section 4. Thus, we return to the discussion of the set of 
equations (4.17) and (4.18). In the light of the above discussion it is justified 
to replace the nonlinear terms in (4.17) and (4.18) (i.e. the first terms on the 
right-hand sides) by their fundamental harmonics. 

A Fourier analysis shows that if uzy is sinusoidal, as defined by (4.15), 
the fundamental harmonic of the nonlinear term 


—im' QU [Us F ur 


: U 
—Ón'0:o'üg'g (=) : cos(o * t), (5.12) 


III 


where 


2 
[2x arsino +5 TR 0e» ifO0<x <1 


fi<x. 


g(x) = 


x ale 


v 


(5.13) 


The graph of the function g is shown in Fig. 5.2. 

Replacing now the first term on the right-hand side of equation (4.17) 
by the expression (5.12) and the first term on the right-hand side of 
equation (4.18) by the corresponding fundamental harmonic leads to the 
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Figure 5.2 f : 
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following simplified equations for the fundamental harmonics: 
Agi e? Apc 


(1+ sp) * Lg o gg + y Ur a 


U 
: cos(o ` [f — tg) 6m 0 tin "8 (32) ` cos(o > t) 
Unt 


— Sm ` Lim’ © ` Un, (5.14) 
deep Amc? 
V Ur V Up 


dts): Lr tip + 


U 
=% o ú? (Z2) costo oft =a) -sr Ly ot (5.15) 
I 


A procedure to solve the system of equations (5.15) and (5.16) is given in 
the Appendix. 

The lengths of the tube elements (particularly Lzy) used in the relations 
(A3) to (A6) should not correspond to the geometric lengths, but to the 
geometric lengths including a length correction due to the particular 
impedances of the tube ends (see Morse and Ingard [10], section 9.1). For 
an unflanged open end, for example, the length correction corresponds to 
about 30% of the tube diameter (see Levine and Schwinger [11]). A further 
virtual increase of the tube lengths is due to higher orders of wL,/c, oL jc 
and ®L7,,/c in (Al) which are ignored within the Helmholtz approximation. 
An analysis which includes these nonlinear terms leads to the “effective 
length" 


1 /o\ 
Lig, Lm: [ T 3 (2) Lc Lin) | (5.16) 
for the situation illustrated by Fig. 5.1, where Lzy is the length which 
already includes the “impedance correction" and L;;,,, also includes the 
first correction to the Helmholtz approximation. 
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6. The solutions and their properties 


Increasing the mass flow through a Helmholtz attenuator as shown in 
Fig. 1.1 (and keeping the amplitude of the pressure oscillation in the 
pulsating system constant) leads to a reduction of the velocity amplitude ù; 
and consequently also to a reduction of the sound attenuation. Further- 
more, the resonance frequency is reduced and the width of the resonance 
band increased. 

Table 6.1 defines a set of data which is used as a case of reference. 

Note that for this reference case, which corresponds to an experimental 
situation discussed subsequently, there is no real pulsating system down- 
stream of the attenuator, but a tube of cross-sectional area A, which is 
also air-filled and the oscillation is generated by a loudspeaker (see next 
section). 

The ratio rp, of pressure amplitudes (i.e. pressure amplitude in the 
resonator vessel divided by f) is obtained with the help of a simple 
calculation, making use of the relations (4.14—(4.16), 

2 
Fpa = ss. (Am ` ty + Az: ür cos(o + o]? 
o'V- P 


+ [4; dr sin(@ : 1) 2. (6.1) 
The mean attenuation power of the Helmholtz resonator is 


Wa =z: Arm’ ûn’ P coslo - to), (6.2) 


and the power of a running wave in a cylindric tube of cross-sectional area 


Ar, Which produces a pressure amplitude 5 upon reflection at a rigid 


Table 6.1 

Data set defining the reference case 

pressure amplitude p — 2.1 mbar 
resonance frequency according to equation 2.13 w —-2:n-: 217.5 Hz 
mean velocity in attenuator tube U = 0; 4 m/s 
mean pressure (pr) = 0.97 bar 
air temperature T=293K 

loss coefficient of inlet tube (2.5 

length of inlet tube L, = 100 mm 
diameter of inlet tube D, = 2.8 mm 
volume of attenuator vessel V — 0.0381 

loss coefficient of attenuator tube Cry = 0.4 

length of attenuator tube £y,100 mm 
diameter of attenuator tube Dy, = 9mm 
length of pulsation tube Ly = 1200 mm 


cross-sectional area of pulsation tube 
gas temperature in pulsation tube 


A, = 0.005027 m? 
T,-298K 
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endwall, is 
Ar: p f 


W,- =, 
ý 8:Qr'Cr 


(6.3) 
where gr and c; denote density and sound speed of the gas in the pulsation 
tube. Thus, it is assumed that the amplitude of the pressure oscillation 
exhibits a maximum at the rigid endwall. 

The relative attenuation of the Helmholtz resonator is defined to be 


r Wa 
A= IL 
W; 


Figures 6-1—6.4 show rp4, r4, yy and (o * fo) as functions of the frequency 
o /2x for the data set defined by Table 6.1. 

A comparison of the results for the cases with and without flushing 
confirms the qualitative statements at the beginning of this section. 

To illustrate the influence of the thermoacoustic boundary layers in the 
tube elements of the resonator 5 and (Pr) are increased by the factor 10. As 
a consequence v becomes 10 times smaller and s; and szy are decreased by 
the factor 10. Comparing the results in Figs. 6.1 and 6.2 shows that for 
weaker boundary layer effects the ratio of amplitudes and the relative 
attenuation are increased. Furthermore, the plot of the phase angle œ ` t; 
shows (see Fig. 6.4) that the resonance frequency is also increased. 

The loss coefficients ¢; and yy deserve our special attention. Increasing 
Cm by the factor 2 with respect to the case having Uj; —4 m/s, (Pr) = 
0.97 bar leads to a reduction of rp, of about 20% and a reduction of the 


(6.4) 


— U n70 «p,» =0.97bar, p^ = 2.1mbar 
-—- Uy =4m/s, «p prO9ten p= 2.1mbar 
----- U = 4mis, Pp 9.7bar, — p^ = 2imbar 


< 20 
„a 

© 15 

oO 

= 

& 

& 10 

4 

$6 5 

2 

E 

E. o 

180 200 220 240 260 


Frequency [Hz] 


Figure 6.1 
Ratio of amplitudes as a function of frequency at different flushing velocities and pressure levels. 
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Figure 6.2 
Relative attenuation as a function of frequency at different flushing velocities and pressure levels. 
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Figure 6.3 
Amplitude of velocity oscillation in the attenuation tube as a function of frequency at different flushing 
velocities and pressure levels. 


resonance frequency of a little more than 1%. Dividing 65, by 2 has the 
opposite effect. 

Increasing C; and keeping the pressure drop across the resonator and the 
mass flow constant, which implies that A, has to be increased correspond- 
ingly, leads to a strong reduction of rp, and also to a strong increase of the 
resonance frequency. 
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Figure 6.4 
Phase angle as a function of frequency at different flushing velocities and pressure levels. 


Hence, it is quite important to keep both loss coefficients as small as 
possible. The values of (;; are particularly critical. By choosing rounded 
tube ends (see Fig. 1.1) flow separation in the tube element JII can be 
largely suppressed (see Sturtevant and Keller [12]), provided the following 
Strouhal number criterion is satisfied 


[() Ry 1 
ee OS € 6.5 
(s. Ta) = ( ) 


where Ry is the radius of curvature of the rounded tube end. 

A last aspect which deserves a few additional remarks is associated with 
the nonlinearity of Helmholtz resonators. As a consequence of the nonlin- 
earity (see e.g. equations (5.14) and (5.15)) 


Um ~ pl? (6.6) 


at resonance, provided the effects of the boundary layers remain weak. 
Hence, according to equation (6.2) 


Wa p". (6.7) 
On the other hand, according to equation (5.3), 

W; ~ p’. (6.8) 
Thus, we obtain from (6.4), making use of (6.7) and (6.8), 


1 
pee, (6.9) 
A pi? 
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which implies that the relative attenuation at resonance becomes weaker 
with increasing amplitude of the pressure oscillation in the pulsation tube. 
For this reason, the attenuators should be designed with respect to the 
largest amplitudes of the pressure oscillations expected in the pulsation tube 
or in a combustion facility. 


7. Experiments 


The goals of the experimental investigation are the verification of the 
analytical model, the acquisition of values for a priori unknown parameters 
like the loss coefficients and the demonstration of the sensitivity on flushing 
and on the shape of the throat ends. This was done in a model attenuator 
using harmonic single frequency excitation. 


7.1. Experimental setup 


The pulsations are generated in a tube which is equipped with a 
loudspeaker on one end and a Helmholtz resonator on the opposite end 
(Fig. 7.1). The dimensions are given in Table 6.1. Large pulsation ampli- 
tudes can be achieved if the resonance frequencies of the tube and the 
resonator coincide. For an arrangement with an open and a closed end the 
tube length has to be 4/4, 34/4, 52/4, ... . The resonator is connected to an 
air plenum which supplies the flushing air. The air flow is determined by 
measuring the pressure difference between plenum and tube as well as the 
temperature in the plenum. 

Microphones are used to detect the pulsations in the tube and in the 
resonator. They are flush mounted in the tube wall and they are embedded 
in absorbing materials (e.g. foamed material) in order to avoid solid body 
sound transmission. 


microphones 


loud 
speaker d resonator 


Figure 7.1 
Helmholtz attenuator test rig. 
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The performance of the resonator is determined by two different meth- 
ods. For the first one the pulsations are measured at two of the positions x,, 
X5 or x; in the tube (Fig. 7.1). For description we choose the positions x, 
and x;: 


1 
Api = A el) = p al fm) — g(x); (7.1) 


Ap, = Be'o- 99-7 of fxs) ~ go) (7.2) 


f and g are the upstream and downstream propagating waves: 


fa exp} f of = En) + al (7.3) 
g=b apid o G ATE iix) + 2 (7.4) 


with amplitudes a, b and phase angles a, f. 

Having the amplitudes A, B and the phase difference ; — ®, of the 
pulsations the real part of the reflection coefficient can be calculated which 
is the ratio of the amplitudes of the downstream (reflected) and the 
upstream propagating (arriving) waves g, f: 


BV B A(X — xi) 

; b 2 (2) 27 e (0.0) - "Im rem] 
i -(?) as BV PEZ pog w(x. — xi) M 

(a) rim eme | 


M is the Mach number of the flow in the pulsation tube. In the present 
investigation M was always very low because the axial flow was only due to 
a small amount of flushing air. Since there is no further source of attenua- 
tion in the present arrangement (if we neglect the small contribution of the 
thermoacoustic boundary layers in the pulsation tube), the reflection co- 
efficient is a direct measure for the attenuation power of the Helmholtz 
resonator. The attenuation power which is already defined by equation (6.2) 
can be written in an alternative form in the limit M «1, provided the 
Helmholtz resonator is connected to the pulsation tube at the location of 
maximum amplitude of the pressure oscillation, 


ATOTCT 2 2 Arp’ 1 Ex 
W, = =- ——— 7.6 
A 8 (a b ) 20507. l "n r’ ( ) 
where 
P _at+tb_ Ir (7.7) 
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Using equations (7.6) and (6.3) we obtain a relation between the relative 
attenuation (equation (6.4)) and the reflection coefficient: 
NR Led 

Wr l4+r 
Measuring the pulsations in the resonator and in front of it (microphones 
#4, #5 in Fig. 7.1) is the second method to check the performance of the 
resonator. The amplitude ratio (equation (6.1)) and the phase difference can 
be directly compared with the prediction of the analytical model. Further- 
more, the attenuation power can be calculated by applying equations 
(6.2—6.4). This can be used to prove the consistency of both measurements. 

As will be shown during the discussion of the experimental results the 
relative attenuation obtained from both methods is different. The measure- 
ments taken from the microphones in the pulsation tube include a small 
contribution due to the thermoacoustic boundary layer and therefore lead 
to a higher attenuation. However, the differences are small. This is con- 
firmed by measurements of the reflection coefficient in the closed tube 
without attenuator. The values are always very close to one. 

Two techniques have been used to measure amplitude ratio and phase 
difference of harmonic signals. The first method uses an oscilloscope where 
the signals are displayed in x-y mode. Whenever the phase difference 
deviates from 0 or m, which is an indication of attenuation, the signals 
trace an ellipse on the screen. Evaluating the amplitudes and the intersection 
lengths of the ellipse with the x, y-axis allows to calculate the phase 
difference. This method is easy to handle, accurate and provides an impres- 
sive visualization of the attenuation effect. Difficulties occur for not purely 
harmonic or noisy signals. 

The second method makes use of a two channel spectral analyzer. The 
frequency response measurement gives directly the data needed. Because of 
an averaging option it can also be applied to noisy signals. Good agreement 
has been found when comparing both methods. 


(7.8) 


7.2. Experimental results 


In a number of experiments the influence of different parameters has 
been investigated. The parameter variations are summarized in Table 7.1. 


Table 7.1 
Summary of parameter variations 


Model attenuator 


radius of curvature Ry, [mm] 0; 2.25; 15 
diameter of air inlet tube D, [mm] 0; 2.8 
length of air inlet tube L, [mm] 2; 100 


mean velocity in the throat U,,,{m/s] 0;4 
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Table 7.2 contains the full set of experiments. It shows the geometrical 
data and the main experimental and analytical (computational) results. The 
computations are fitted to the experiments with respect to resonance fre- 
quency, amplitude ratio and/or relative attenuation by choosing appropriate 
loss coefficients and corrections of the throat length. fọ listed in Table 7.2 is 
the resonance frequency according to the Helmholtz formula (2.13), 
whereby the Levine-Schwinger correction is taken into account to define 
Lj. For resonators having one throat only also the correction due to higher 
orders (equation 5.16) is applied. 

We first discuss experimental and analytical results for the case of a 
closed resonator (1.e. no inlet for flushing air) with sharp edged throat ends. 
Figures 7.2 and 7.3 show the amplitude ratio rp, (pressure amplitude in the 
resonator divided by the pressure amplitude in the tube, cf. equation (6.1)), 
the phase angle between the two pressure signals and the relative attenua- 
tion r, for experiment No. 1 (El). The amplitude ratio and the phase angle 
can almost perfectly be predicted by choosing a loss coefficient of “yy = 1.85 
and a throat length of Lj; -— 100mm which corresponds exactly to the 
geometrical length. For a sharp edged throat we expect an almost total loss 
of kinetic energy at outflow and an additional loss at inflow into the throat 
due to flow separation (i.e. “jet contraction"). Hence the value of the loss 
coefficient seems to be reasonable. 

In a closed resonator the pulsation in the resonator volume and the 
velocity in the throat exhibit a phase shift of 90? (cf. equation (4.14)). 
Hence, the phase angle between the pulsation in the tube at the throat exit 
plane and in the resonator volume is related to ct; by (cf. equations (4.10) 


Ámpl. rátio: 
--@-- Exper, | 
Theor. 


PA 


Amplitude ratio r 
Phase angle pop [°] 
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Figure 7.2 
Amplitude ratio and phase angle of the closed resonator with sharp edged throat ends (experiment 
E1). 
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Figure 7.3 
Relative attenuation of the closed resonator with sharp edged throat ends (experiment EI). 


and (4.15)): 
o, = 6, = wto a 90°. (7.9) 


At the resonance frequency, f, = 209.5 Hz, the amplitude ratio reaches a 
maximum of about 11 and the phase shift ®&,— ©, crosses —90°, i.e. 
cos(@fo) = 1. In other words, the velocity in the throat and the driving 
pulsation in the tube are in phase which is a necessary condition for 
maximum attenuation (cf. equation (6.2)). At low frequencies, f/f, < 1, the 
pressure in the tube and in the resonator oscillate in phase while at high 
frequencies, f/f, > 1, these pressures oscillate in counter-phase. Due to a 
small effect of the boundary layers the values measured for attenuation in 
the tube are always slightly higher than the computed ones (Fig. 7.3). 
The next set of experiments illustrates the influence of flushing on the 
performance of a Helmholtz resonator which is important for some applica- 
tions. The relative attenuation of a closed resonator, a resonator with air 
inlet tube but no flushing and a flushed resonator (E3, E4, E5) is compared 
in Fig. 7.4. The increase in resonance frequency of about 10 Hz in E4 
compared to E3 agrees very well with the Helmholtz formula for a res- 
onator with two throats (equation (2.13)). The peak attenuation is reduced 
by about 10%. The reason can be seen from equation (A4) by considering 
the resonant case where cos(wf)) = 1. If a second throat is added the second 
term on the right-hand side becomes negative because sin(wt) < 0 as defined 
by equation (A6). A decreasing right-hand side of equation (A4) leads to a 
reduction of yr. As a consequence also the attenuation is reduced. Only for 
an absolutely loss-free air inlet tube (s; = 0, £, = 0) the performance of the 
resonator would not deteriorate. In this case sin(wt) = 0 would be obtained 
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Figure 7.4 
Influence of flushing on the relative attenuation of the resonator with sharp edged throat ends 
(experiments E3, E4, E5). 


from equation (A6), i.e. the gas columns in the two throats would be 
oscillating in counter-phase like the pistons in a boxer engine. 

A further reduction in performance occurs if the resonator is flushed as 
in E5 where U;, —4 m/s. This behaviour appears as a result of the en- 
hanced losses in the throats due to the change of the velocity level, which is 
accounted for by the function g (equation (5.13), Fig. 5.2). Equation (A4) 
is directly affected on the left-hand side and indirectly on the right via 
sin(wt), which assumes a larger negative value. The values which are shown 
in Table 7.2 indicate that flushing shifts sin(wt) towards — 1. The pressure 
in the resonator and the velocity in the inlet tube are then oscillating in 
phase, i.e. the inlet tube acts like a leak in the resonator. 

The experimental results of E3 and E4 can be predicted with only minor 
changes of £j; and Lj; (Table 7.2). In the case of flushing, however, Cry has 
to be tripled and Lzy has to be significantly reduced. The very low Reynolds 
number of the throat flow and asymmetric flow conditions may serve as an 
explanation. Errors in the flow rate measurement cannot explain these 
changes, because Uy would have to be increased by the factor 3.5 to reduce 
the attenuation to the measured value, if the value of £5, is kept to be 1.85. 

Similar experiments as before have been performed for a throat with 
rounded ends (E6, E7 and E8 shown in Fig. 7.5). An arc of 15 mm radius 
and an angle of 40? has been used (Fig. 7.6). In this case the peak 
attenuation increases by more than 55% compared to the sharp-edged 
throat. The resonance frequency undergoes a shift of 6-7 Hz. The higher 
values are due to a reduction of the effective length of about 6 mm. The loss 
coefficients £5, are decreased to 25-40% of the values for sharp edges. 
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Influence of flushing on the relative attenuation of the resonator with rounded throat ends (experiments 
E6, E7, E8). 


Figure 7.6 
Resonator throat with rounded ends used in the experiments E6—ES8. 


Again the flushed case has an increased loss coefficient although the 
Strouhal number is about 0.4. 

A further experiment (E2) was carried out with a moderate arc radius of 
2.25 mm which corresponds to a quarter of the throat diameter. Compared 
to sharp edges (El, E3) the loss coefficient is reduced by 40% and the 
attenuation is increased by 20%. 

The experiments E9 and E10 (Table 7.2) show the effect of replacing the 
air inlet tube (L; = 100 mm) by a simple hole in the resonator (L; = 2 mm). 
The gas column in the hole has practically no inertia such that sin(wr) is 
very close to —1. As already pointed out, this is the situation of a leaking 
resonator. The attenuation is almost completely lost. In this case flushing 
has some "sealing" effect on the resonator because it reduces the amplitude 
of the velocity oscillation 4; in the hole (cf. equation (A6)). 
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Figure 7.7 
Loss coefficient as a function of Strouhal number and flushing velocity (experiments El-E8). 
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Figure 7.8 


Correction of the throat length as a function of rounding radius and flushing velocity (experiments 
EI—E8). 


As already mentioned appropriate loss coefficients and corrections for 
the throat length have to be chosen in order to fit computational and 
experimental results (Table 7.2). The values are shown in Figs. 7.7 and 7.8 
where the loss coefficient £7; is displayed versus Strouhal number (equation 
(6.5)) and the length correction AL, is displayed versus radius of curvature 
Ry, of the rounded tube ends. Parameter is the flushing velocity Uy. Figure 
7.7 indicates a strongly reduced loss coefficient even for Strouhal numbers 
much lower than recommended by the criterion given in equation (6.5). 
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The throat length (Fig. 7.8) needs no additional correction in case of the 
closed resonator with sharp edges where the effective length according to 
equation (5.16) is applied. Rounded tube ends lead to a reduction of the 
effective throat length as expected. The length correction is about 0.4 times 
the radius of curvature. For a resonator with two throats equation (5.16) is 
not valid anymore. As a consequence, the effective length has to be 
significantly increased in the case without flushing. Flushing, however, 
reduces the correction which is necessary to match computation and exper- 
iment. This points towards strong separation effects because of the asym- 
metric flow conditions and the reduced unsteadiness of the flow expressed 
by Ug; [üg, > 0. The behaviour coincides with the strong increase of the loss 
coefficient. 


8. Conclusions 


(i) The accurate prediction of the frequency-dependent attenuation prop- 
erties of Helmholtz attenuators requires nonlinear effects due to flow 
losses, effects due to thermoacoustic boundary layers and the first 
frequency correction to the “Helmholtz approximation” to be taken 
into account. 

(ii) To achieve high attenuation performance the flow losses have to be 
reduced by rounding the ends of the throat and the air inlet tube. So 
far the loss coefficients and the correction of the throat length have to 
be determined experimentally. 

(ii) In case of a flushed attenuator the loss of performance can be kept 
small if the air inlet tube is not too short and if it has a small cross 
section. A length similar to that of the throat is recommendable. 
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Procedure to solve the system of equations (5.15) and (5.16) 


Inserting the fundamental harmonics (4.15) and (4.16) into the equa- 
tions (5.14) and (5.15) and equating terms of equal phase leads to four 
equations for the four unknowns üy7, ty, à; and Tt: 
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Knowing the pressure drop across the resonator (see Fig. 1.1) we can 
determine the mean mass flow and the mean velocities U, and Uj. 
Eliminating first cos(o - t) from (A3) with the help of (A5) and 
sin(o : T) from (A4) with the help of (A6) and, secondly, eliminating (c - to) 
by combining (A3) and (A4) leads to an equation of the form 


für, dj; parameters) = 0. (A7) 


Furthermore, eliminating (w : t) from (A5) and (A6) we obtain an equation 
of the form 


g(ü,; parameters) — 47; = 0. (A8) 


Finally, eliminating 2j; in (A7) by means of (A8) leads to a single equation 
for à; which can be solved numerically by looking for the smallest positive 
zero. Having found Z, we obtain 45; with the help of (A8), (œ > t) with the 
help of (A5) or (A6) and (o : tọ) with the help of (A3) or (A4). 
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Summary 


In the present paper a nonlinear acoustic theory is proposed, to accurately describe the properties 
of a generalized type of Helmholtz resonators. The theory may be used as a layout tool to design sound 
attenuators for combustion facilities. The present investigation shows that, in addition to the nonlinear 
effects, the effects due to thermoacoustic boundary layers and the first frequency correction to the 
“Helmholtz approximation" should be taken into account, in order to predict the frequency-dependent 
attenuation properties of Helmholtz attenuators with the accuracy required by typical technical 
applications. A series of experiments is presented which is used to validate the theoretical predictions. 


Zusammenfassung 


In der vorliegenden Abhandlung wird eine nichtlineare akustische Theorie vorgeschlagen, mit der 
das Verhalten einer verallgemeinerten Art von Helmholtzresonatoren beschrieben werden kann. Die 
Theorie kann als Werkzeug zur Auslegung von Schalldämpfern für Verbrennungsanlagen benutzt 
werden. Die vorliegende Arbeit zeigt, dass neben den nichtlinearen Effekten auch die Einflüsse der 
thermoakustischen Grenzschichten und der ersten Frequenzkorrektur zur “Helmholtzapproximation” 
berücksichtigt werden müssen, damit die frequenzabhängigen Eigenschaften von Helmholtzdämpfern mit 
technisch ausreichender Genauigkeit beschrieben werden kónnen. Mit einer Reihe von Experimenten 
werden die theoretischen Ergebnisse untermauert. 
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